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Abstract 

We present a simple model for chiral molecules which yields the frequency-dependent 
multipole-multipole polarizabilities required for calculating the dispersive van der 
Waals constants in a simple, closed, and consistent form. The model consists of a 
single effective charge confined to a one-dimensional wire formed by two circular 
arcs, which are twisted with respect to each other by an angle characterizing the 
chirality. The resulting polarizabilities show a simple functional dependence on the 
model parameters, and they serve to mimic the chiral properties of small molecules, 
such as H2S2, in a reasonably realistic way. 
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1 Introduction 

Wire models [T1I2] are popular to describe molecular configurations where elec- 
trons are allowed to move almost freely along quasi one- dimensional structures, 
as is the case e.g. with the dislocated 7r-electrons in organic chain molecules 
[3f4f5] . Studies comparing these models with ab-initio electronic configuration 
calculations [2|6] show that, in spite of their simplicity, these models succeed 
in reproducing many qualitative features and yield even quantitatively rea- 
sonable results. 

In the present article we discuss a simple wire model which aims at describing 
small chiral molecules, and in particular the dependence of the chiral nature 
on their electronic susceptibilities. The latter determine important molecular 
properties such as the optical activity and the dispersion forces acting between 
molecules. While the chiroptical properties have already been subject of several 
publications [32EI1E], the present article is focused on the dispersion forces. 
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In general, it is a difficult task to calculate the dispersive interaction poten- 
tial strengths between chiral molecules (with few or no symmetries) and other 
polarizable particles if one has to go beyond the bulk van-der-Waals force 
U = —C 6 /r 6 . This is the case, e.g. if one needs to assess the relative strength 
of chirality- discriminative part of the interaction versus the non-discriminative 
ones. The reason is that multipole-multipole polarizabilities are required on 
the whole imaginary frequency spectrum for calculating the higher order dis- 
persion constants which are linear combinations of integrals [7] 

where aff'^L/ and aff'^f,/ denote the polarizabilities for molecules A and B, 
respectively, in spherical tensor notation, see Eq. (J5TJ) below. 

The van-der-Waals dispersion constant Cq depends on the electric-dipole- 
electric-dipole polarizabilities of both molecules {I a = £'a — = ^'b = !)■ 
While it does not discriminate between the left- and right-handed configu- 
rations of an oriented chiral molecule, higher order dispersion interactions 
differ in general for left and right configurations. Often, a rotational average 
over all orientations of the chiral molecule cancels this left-right distinction in 
the cross sections. However, consider a coherent superposition of the left- and 
right-handed configuration states [8], which will be subject to decoherence due 
to collisions with the atoms of a surrounding gas. In this case, the decoherence 
process is determined by collisions of spherical gas atoms B off the oriented 
chiral molecule A. Here, the lowest order discriminative interaction with B 
(£b — £'b = 1) ^ s gi ven by the electric-dipole-electric-quadrupole polarizability 
[9J, £a — 1 and £' A = 2. The resulting dispersive interaction depends on the 
distance r between the molecules as r~ 7 [10]. We will focus on this dominant 
discriminative interaction below, motivated by the above-mentioned question 
of collisional decoherence. We note that there might be additional forces of 
comparable strength which do not discriminate left and right, and that only 
non-retarded forces are considered in the following. For a discussion of retarded 
interaction between chiral molecules, see [TTlll2|li3| . 

The main aim of the present model is therefore to provide a minimalist, while 
consistent and trustworthy model for the bulk electronic properties of chi- 
ral molecules. In particular, we give explicit and closed results for the rota- 
tory power, the electric-dipole-electric-dipole polarizability, and the electric- 
dipole-electric-quadrupole polarizability, both at imaginary frequencies. The 
simplicity of the model even admits the calculation of further higher multipole- 
multipole polarizabilities with small effort. 

The model is specified by two parameters, the length L determining the bulk 
polarizability and the twist angle Xi which corresponds to the dihedral angle 
in a molecule, see Fig. [TJ By choosing the parameters to represent the chiral di- 
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hydrogen disulfide, H2S2, we obtain reasonable values for its various electronic 
properties, suggesting that the model may serve as a good approximation for 
a typical small chiral molecule. Yet, the aim here is of course not to reproduce 
precise numbers, but to have an easily accessible, consistent and plausible de- 
scription of the dependence of the electronic susceptibilities on the molecular 
parameters, and in particular on the dihedral angle. 

The structure of the article is as follows. Sect. [2] presents the model and its 
quantization based on the canonical description of the one-dimensional dynam- 
ics. The explicit forms of the most important Cartesian operators (position, 
momentum, magnetic dipole, and electric quadrupole moment) are then de- 
rived in Sect. [31 together with closed expressions for their matrix elements in 
the energy eigenbasis. Based on these results the rotatory strength is calculated 
in Sect. SJ as a function of the chiral angle. Sect. contains the calculation 
and discussion of the electric-dipole-electric-dipole and electric-dipole-electric- 
quadrupole polarizability for imaginary frequencies. An example for interac- 
tion strengths derived from the polarizabilities is discussed in Sect. [6] and we 
present our conclusions in Sect. [3 




Fig. 1. The twisted arc model is given by a wire of length L in the form of two arcs 
formed by quarter circles. They can be twisted with respect to each other by an 
angle %■ This way a three-dimensional chiral structure is formed, except for \ = 
and x = 7r, when the wire is in its planar trans- and its cis-configuration. 
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2 The twisted arc model 



Our model aims at describing the excitation properties of a chiral quantum 
system in the easiest possible way. It is formed by two connected circular arcs, 
C\ and Cu, each described by the common radius R and an angle of 90°. The 
total length of the wire is thus L = irR. We put the origin of the coordinate 
system at the junction of the arcs such that C\ lies in the xy-plane. For \ — 
also Cu lies in the xy-p\a,ne (trans-configuration), while for x 7^ it is turned 
around the x-axis, see Fig. |5J We call x the twist angle, and take it as the 
angle between the osculating plane of Cu and the y-axis (— n < x ^ For 
X = 7r the wire is in its cis-configuration, and for x ^ 0, it the angles x an d — X 
correspond to configurations with opposite chirality. Alternatively one could 
choose the dihedral angle <p = tt — x- 



Fig. 2. Projections of the twisted wire model onto the rry-plane (left) and the 
yz-plane (right). The arcs C\ and Cu are quarter circles with a radius of R = L/ir. 
C\ lies in the xy-plane as indicated and the twist angle x is defined to be positive 
(negative) for Cu pointing into the positive (negative) ^-direction. The origin of the 
parameterization coordinate s is taken to be the loose end of C\. The C2 symmetry 
axis, drawn as a dashed line, lies in the yz-pane. 

A charged particle of mass M and charge q is confined to move freely along 
the wire defined by the two arcs. It is constrained by an infinite potential step 
at both ends, but there is no force in the x-direction at the joint of C\ and Cu- 

The model shows C2 symmetry. The symmetry axis (dashed line in Fig. [2]) lies 
in the wz-plane and includes the origin. It has an angle of — x/2 with respect 
to the z-axis. 

We note that an alternative construction, suggested in [5] , would be to take (at 
least) three straight wire parts joint together at angles of 90°. We prefer the 
present cornerless structure, since it consists only of two equivalent legs which 
mirror the symmetry of the wave function, and admit simple closed expressions 
for the relevant matrix elements, see below. For a critical examination of the 
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thin wire model with corners see [13]. Another widely used model for chiral 
molecules is the helical thin wire model [15], which however does not allow 
the straightforward modelling of the dihedral angle of H2S2. 



2. 1 Quantization of the twisted arc 



In order to quantize the motion let us first consider the canonical description 
of the classical motion p] . The natural generalized coordinate is the length of 
the wire s (0 ^ s ^ L). Noting R = L/n we find from Fig. [2] its relation to 
the Cartesian coordinates. 
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As can be easily checked, this parameterization guarantees that the kinetic 
energy T is given by the generalized velocity 



T v x +y + z ) = Y s ' ^ 

In the absence of a (magnetic) potential the Lagrange function equals the 
kinetic energy, £ = T. The conjugate momentum is defined as p s = dC/ds = 
Ms, so that the Hamilton function H = T takes the form H = p 2 s / (2M). The 
Cartesian components of the momentum are thus obtained, via p x = Mx = 
M(dx/ds)s, as 
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We turn to the quantum description by replacing p s by the differential operator 

acting in L 2 ([0,L]). Operators will be indicated by sans-serifs throughout. 
The Hamiltonian is now 

h 2 d 2 

H - ~w <«) 

Together with the boundary condition ip(0) = ip(L) = 0, it yields the eigen- 
functions of a particle in a one dimensional box with length L 

2 /mrs\ 

ipn{s) = y-^sm I — 1 ,n = 1,2,3, . . .. (7) 

The corresponding energies depend quadratically on the excitation number, 
E n = n 2 Ei, with the ground state energy given, for M = m e , by 



„ _ 7T 2 ( Ly 2 35.06 eV 

El ~^\VJ haxtree ~ ~{L/kf' (8) 

Choosing L as the sum of the binding lengths of H2S2, i.e. L = 8.6ao, yields a 
first excitation energy of E 2 — Ex = 3Ei = 5.4eV, which is in good agreement 
with the result of ab initio calculations for the excitation energy, 5.48eV, for 
the longest wavelength transition [T6] . 

In the following, we will use atomic units (au), h = eo = m e = 1, with M = m e . 



3 The Cartesian operators 



The functional dependence of the Cartesian operators on the coordinate s will 
in general be different on C\ and C\\. It is therefore convenient to introduce for 
any operator A the corresponding partial operators A^ and A^ 1 ^ satisfying 

A {a) = e (| - a ) A® (s) + 9 (3 - I) A(") (a) . (9) 

The matrix elements of A with respect to the energy eigenstates (I7|) can then 
be calculated as 



(m|A|n) = [ L/2 dsiP* m (s)AV ( s ) ^ n ( s ) (10) 







+ f L dsr m (s)^ n) (s)Ms). 

JL/2 
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The C2 symmetry of the wire model will allow to express matrix elements of 
operator A by the restricted operator A® alone, see below. 



3. 1 The position operator 

The Cartesian dipole operator is determined by the Cartesian components of 
the position operator, d = q~x . From ([2]) we obtain immediately 
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As an advantage of the present model, the matrix elements can be shown to 
be simply interrelated, 



(m|z (II) |n) = (-) 
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It follows that only a small fraction of the matrix elements needs to be evalu- 
ated. For those we obtain 



(m |x| n) = (m nj + (m 
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Here we introduced abbreviations which will be used throughout the paper, 
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2 ^,2 



D mn :=m 2 -n 2 . (14) 



Moreover, the following factors show up frequently: 
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(15) 



(16) 



The matrix elements of dipole operators in spherical form are calculated from 
the basic matrix elements Eq. ( jl~3l) 
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TTie momentum operator 



The Cartesian components of the momentum operator on the individual arcs 
can be obtained for the corresponding classical expressions (J3J). However, un- 
like the full operator p, the hermiticity of and is not guaranteed. 



By replacing the canonical momentum p s in the symmetrized version of (jlj) 
by its corresponding operator dSJ), and carrying out the derivatives as far as 
possible, we get the Cartesian momentum operators 
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The operators p^, pQ 1 ', p^\ and are hermitian, while pj|P and p^ 1 ', are 
individually not hermitian, but only their sum. For example, p)|P obeys the 
relation 

* 2 

n|p®|m) = ^m|p®|n) ~ ^isnSm- (19) 

The second term on the right hand side is the boundary term of the partial 
integration showing up if the differential operator is applied to the bra. It 
does not necessarily vanish here, since the wave function may be finite at the 
boundary point s = L/2. A remedy would be to define a modified momentum 
operator 

- ^) ■ (20) 
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This operator on C\ is now hermitian, (m\p^\n) = (ra|pW|m)*, since the 5- 
function cancels the contribution of the boundary term. Similarly, the defini- 
tion 

PP := Pi H) " (s - (21) 

ensures both the hermiticity on C\\ and the relation p x = pW + pP 1 ). In general, 
it would be therefore more convenient to express momentum matrix elements 
in terms of the hermitian operators (T2"U1) and f[2"Tj) . However, below the momen- 
tum operator will occur only as part of the magnetic dipole operator, where 
this modification is not required, as discussed next. 



3.3 Magnetic dipole moments 

The magnetic moment is proportional to the angular momentum operator 
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where we still use atomic units and allow for an effective charge q = Z e geo, 
which might be useful fitting parameter. 



The magnetic dipole moment depends on the origin, and it is natural to choose 
the position of the joint between the two arcs in the present model. It can be 
easily seen that the non-hermiticity of the parts and p^ is then irrelevant 
since the x-component of the momentum is multiplied by y or z, which vanish 
at s = L/2 and thus suppress the hermiticity terms in Eqs. (1201 and ff2T]) . The 
ability to evaluate the magnetic moment directly, without having to resort to 
hermitian corrections, is another advantage of the present model. 

Again one can derive helpful inter-relations of the matrix elements, 



m\ (xp y - yp x ) (II) |?^ = (-) n+m cosx (m\ {xp y - yp x ) (1) \nj (23) 
{m\zp x - xp z \n) = (m\ (zp x - xp 2 ) (II) \n) 



) smx\m\ [xpy -yp x ) KJ \n 
(ii) 



(m\yp z - zp y \n) = (m\ (yp z - zp y ) ( \nj = 0. 

The last relation follows immediately from the proportionality of y and z 
components of the operators and the vanishing of the z components on C\. It 
follows that a single integral remains to be calculated. 
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to those defined in ([T4"]) - (IT6"]) . The hermiticity of the part is now evident, 
(n|m®|m) = (m|m^|n)*, implying with (123]) that all the components of rrt 
are hermitian. 
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3.4 Electric quadrupole moments 



The quadrupole moments are given by a quadratic combination of the position 
operators. In Cartesian coordinates, r = (x,y,z), we have 



(27) 



Again, a variety of proportionality relations serves to reduce the calculational 
effort considerably, 
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It follows that three integrals need to be evaluated, 
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The last two expressions, (!30|) and (13TT) . are valid for those combinations of 
values n and m for which the denominators do not vanish. We omit the general 
results for m—n = 1,2 since the matrix elements will be needed only for m = 1 
below. In this case they take the form 




Using these results one gets quite compact expressions for the quadrupole mo- 
ments, which display a simple dependence on the twist angle \. Here we note 
the matrix elements of the spherical quadrupole operators, Q2,^, with magnetic 
quantum number /x. They are required for the calculation of the electric-dipole- 
electric-quadrupole polarizability in tensorial form, a k ' k , (u>), see below 




12 



4 The rotatory strength 



We proceed with the discussion of the rotatory strength, the most prominent 
property of chiral molecules. We do so to demonstrate the ability of the twisted 
arc model to display electromagnetic properties of real, chiral molecules. The 
latter are usually optically active, and their chiroptical properties depend on 
the optical rotatory tensor |17j . 

i^Im^lldilnKnKIl)}. (34) 
The trace of the optical rotatory tensor yields the rotatory strength, 

R nl = Rf x + Rf y + Rf z = Im {(l|dV) ■ (n| m |l}} . (35) 

It determines the optical rotation angle for orientationally averaged molecules. 
Recognizing an energy-resolution of the identity operator one finds that the 
sum over all states must vanish, 

oo oo 

Y,R nl = Im£ (l| d \n) ■ (n|m|l) = Im(l| d • m |l) = 0, (36) 

n=l n=l 

because the expectation value of the hermitian operator d • rn must be real. 
We note also that the rotatory strength does not depend on the choice of 
origin. The most convenient choice of origin is therefore the joint of the arcs, 
since hermiticity corrections of the momentum operator are then not required, 
as discussed above. Note that an object independent of the origin is obtained 
by amending the rotatory tensor (|34|) with additional terms depending on the 
quadrupole-dipole tensor, which however cancel out after rotational averaging 

From the above formulas for the matrix elements of rn and d one obtains 
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for n > 2, 



(38 



As one expects, the rotatory strength is antisymmetric with respect to the 
twist angle, and the sum rule (1361) for the rotatory strength is fulfilled. 



Reasonable parameters for a small chiral molecule such as for H2S2 are L = 8.6 
and Z 2 S = 4.5 (see Sect. 15.21) . With these values one finds that R 21 for x — f 
is 0.388 au, corresponding to about 180 x 10~ 40 esu, which is larger by a factor 
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of 5 to 10 compared with theoretical values for H2S2 [T6][l~9] . This is mainly 
due to cancellation effects generated by the two lowest lying excited states 
of H 2 S2, which are nearly degenerated for dihedral angles around it/2. They 
have rotatory strengths which are similar in absolute values but opposite in 
sign [T6] . 



Having found reasonable choices for the parameters L and Z c n, we now proceed 
to evaluate the frequency dependent rotatory power G' (u), a further quantity 
which permits a comparison with the literature values of the real molecule 
H 2 S 2 . The rotatory power is expressed by the electric-dipole-magnetic-dipole 
polarizability [20] . 

00 B"! 

G' n (u) = -2lu y ^ . (39) 

IJ V / ^ I 771 771 \2 2 V ' 

n=2 - hi) -U z 

The average over all orientations of G'^ is given by the trace 
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G»= E G« (w) = -2w S; — ^ 2 ~ , (40) 

i=x,y,z n=2 



(En-Eif-U 2 



which can be expanded for small frequencies (compared to the excitation gap). 
For x — 7r /2 we find 

G' w/2 (lo) Z^L- 3 = -1.142 ■ 10- 4 lu - 1.61 • 10~V + O (a/) . (41) 

Here, the frequencies are in units of Ei/h. The polarizability, given in atomic 
units, scales as Z^ S L 3 , where a factor Z^ a L is contributed by the dimensions 
of the rotatory tensor and a factor L 2 by the decrease of the excitation energy 
with increasing length scale. 



The specific rotation angle per dm can be calculated as (see e.g. [TO] ) 

^H.-W.S^-sinx^l^g), (42, 

where rj is the concentration of chiral molecules modeled by twisted wires. The 
validity of ( 1421) is restricted to frequencies well below the first excitation energy. 
Inserting the frequency for the sodium D-line (i?Na ~ 1.17J5i for L = 8.6ao), 
with the above choice of L = 8.6a and Z 2 S = 4.5, one gets 

0Na-D = -103° dm" 1 sin x ^ mol , ^ ■ (43) 

The specific rotation angle of H2S2 has a somewhat different functional depen- 
dency on x due to abovementioned cancellation effects, but its value is about 
the same order of magnitude as in CI calculations [I~9| . 
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Note that the optical rotatory tensor (|34|) appears in the discriminatory part 
of the dispersion interaction between two chiral molecules (here without rota- 
tional averaging) [TT] . 

E ch _ ch = — (fo - 3RA) fe - 3RA) £ p l \ % ■ (44) 

" n.m f-'nl ~r -'-'ml 

We will not evaluate the expression further since we are mainly concerned 
with interactions between a chiral molecule and achiral atoms. 



5 Electric multipole-multipole polarizabilities 



We proceed to calculate the electric multipole-multipole polarizability tensor 
required for evaluating the dispersion interaction coefficients, Cq and C-j. With 
this application in mind, the results will not be presented in full generality, but 
will be confined to the ground state polarizability at imaginary frequencies. 



5.1 General form of the electric- dipole- electric- dipole polarizability tensor 



Apart from a xx , all elements of the polarizability tensor ct depend on the twist 
angle x- We start by discussing this dependence in more detail, noting that 
the elements of the polarization tensor of the state m are given as 
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where i, j G {x, y, z}. 



Due to the time-reversal invariance the tensor elements are real at imaginary 
frequencies, see e.g. [21], implying that the tensor is a symmetric matrix and 
can hence be written as 
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For the special case m = 1 and inserting d = Z e ^T we get finally 
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where we absorbed a factor n 2 /2L 2 into the frequency, thus measuring the 
frequency ui again in units of Ei/h, see (151) . The index indicating the reference 
state m — \ will be suppressed from now on. 



We find that, due to the interrelations of the position matrix elements noted 
above, the frequency dependence of the spherical tensor is determined by only 
three functions, 




With their help the structure of a at imaginary frequencies is readily specified, 



a (uj) 



|-(l + c x )a(a;) 



lt s x a M 



(49) 



£ (1 + c x ) a (u) g+ (w) (l + c\) + (u) c x (g + (u) + g~ (u) c x ) s y 



h s x a (^) 



{g+ (w) + g- (w) c x ) s x 



9+ \U) K 



Here we used the abbreviations s x := sinx and c x := cosx and g± : = 
g c (u) ± g Q (u). The sums in ( l4"8j) converge rapidly as the terms are of order 
O (n -6 ) for large n. 



The dependence on the twist angle is particularly simple in f)49p due to the 
proportionality relations (Ti~2]) . One observes that only the off-diagonal elements 
in the third row and the third column change their sign when switching to 
opposite handedness, x ~^ ~X- Such a change of sign in x is equivalent to a 
reflection at the xy-plane, z — > —z. After an additional rotation by 7r around 
the z-axis a full parity operation, f — > — r, is obtained. Hence, one gets the 
tensor of the — x configuration after rotating the polarizability tensor of the 
+X configuration, which illustrates the well known fact that all dispersive 
interactions solely derived from the electric-dipole-electric-dipole polarizability 
tensor do not discriminate left- and right-handed enantiomers. 



1(3 



5.2 Choice of parameters 



In the following we present numerical results for a specific choice of parameters, 
which are adapted to mimic the literature values for H2S2, e.g. [T6l[T9] . First we 
fix the length parameter L by the sum of the binding lengths of the molecule 
which is found to be 8.6 do- The twist angle in the wire model is chosen to be 
X = if which is a good approximation of the dihedral angle in the molecule. 

To fix the effective charge Z eff we consider the static (i.e., u = 0) values 
for the electric-dipole-electric-dipole polarizability g + (0) = 1.963 x 1CT 4 au, 
(0) = 1.218 x 10" 4 au, and a (0) = 3.42 x 1CT 3 au. The values of g + and 
</_ as function of uj are depicted in Fig. [3j A comparison with the Lorentzian 
Z 2 R L A a (u) from (I48p . which has a width of 3 and a strength of 84.2, shows 
that the latter dominates the polarizability tensor. 

The rotational average of the static polarizability for \ = ±7r/2 is then given 
by 

L 4 Z 2 

L 4 Z 2 cS a (0) = (a xx (0) + a vy (0) + a zz (0)) (50) 

^L 4 Z e 2 ff 1.271 x l(T 3 au. 

Here the contribution of a xx (0) dominates the average static polarizability. A 
comparison with the literature value [IE] for H2S2 of about 31 au results in 
Z c 2 ff = 4.5 (Ztf «2.12). 




0.5 1 1.5 2 2.5 3 

Fig. 3. Frequency dependence of the functions Z 2 s L 4 g + (io) (solid line) and 
Z 2 s L 4 g~ (ll>) (dashed line) [atomic units; L = 8.6; Z 2 S = 4.5]. The frequency is 
given in terms of the ground state energy E\ / h, see © . 
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5.3 Electric- dipole- electric- quadrupole polarizability 



The dominant discriminative dispersion interaction between an oriented chi- 
ral and an achiral atom or molecule depends on the electric-dipole-electric- 
quadrupole polarizability tensor [SJ. In addition to the dipole matrix elements, 
it is determined by the matrix elements quadratic in the coordinate operators, 
(n \ u r j \ m )i discussed in Sect. 13.41 

For the calculation of dispersion constants it is preferable to consider polariz- 
ability tensors for imaginary frequencies in spherical representation 



a^^) = -^(2t+l)(2t> + l) (51) 

x Yl 7^ — " , 2 \. .2 <i|Q/,aI^> HQ^'I 1 ) , 

n=2 {En -Ex) - (lUj) 

with the spherical multipole functions 



Q e , k (f) = Zcsr^^-^Y^ (6, <p) . (52) 

For the electric-dipole-electric-quadrupole polarizability the multipole opera- 
tors have ranks i = 1 and £' = 2, 

°<%k> ( iu} ) = £ (T? En rp ^\ — 2 ^l dfe l n ) ( n IQ2,fc'|l) • (53) 
2tt n=2 (E n -E^ +u 2 

The required matrix elements can be found in (TPTj) and (j33l) . Obviously, the 

1 2 

polarizability a k ' k , transforms under parity operation like a product of three 
coordinates, 1^17, which results in an overall negative sign. Thus, the disper- 
sive interaction derived from it discriminates left-handed and right-handed 
molecules. Finally, it should be noted that the electric-dipole-electric-dipole 
polarizability a 1 ' 1 does not depend on the choice of the origin, while the 
electric-quadrupole-electric-dipole polarizability a 2 ' 1 does [TU] . 



6 Potential strengths 



As an illustration for the use of the polarizability tensors, let us evaluate the 
interaction potential between a helium atom in ground state and the twisted 
arc. Its calculation requires the electric-dipole-electric-dipole polarizability of 
helium. Since helium is spherically symmetric only a single element of the 
spherical polarizability tensor does not vanish, namely 0^0 . For our purposes 
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it suffices to approximate the helium polarizability by 



a o,o M = — J' ( 54 ) 

with g = 2 and wn e ~ 1.33 (atomic units), which is a good approximation for 
frequencies well below the excitation energy [22]. 

After evaluating the integrals (CE]) the dispersion potentials can be calculated 
following the theory in [23] by disregarding retardation effects. They assume 
the form 

^ 6 (f) = -^p and U 7 (f) = -^P- (55) 

with r the distance vector between helium atom and the center of mass, and 
r = f/r. For simplicity, we take the center of mass to lie in the origin. The 
potential strength C 6 (f) is calculated using the spherical dipole-dipole polar- 
izability tensor, a 1 ^,, of the twisted wire model which is readily obtained from 
the Cartesian polarizability (149|) . The calculation of Cj (f) uses the dipole- 
quadrupole tensor ( 1531 . 



(56) 



To compare the interaction potentials we consider the surface of equality 

0,1 ' \C,(r)\- 



Thus, at distances r = £r eq the strengths of the potentials are related by 
\Uq\ = £\U?\. Since the potential strengths are linear combinations of the 
polarizabilities r cq (f) is proportional to L. At the same time, the dependence 
on Z e s cancels out. 

Figure H] shows C-j/Cq for twist angles x — 71 an d X — ±tt/2 at L = 8.6ao- 
Here, the direction of the helium atom, r = (0,<p), is specified with respect 
to the symmetry axis of the twisted wire (the dashed line in Fig. [2]). The 
azimuthal angle (with respect to the x-axis) is chosen to be if = tt/2. Negative 
values of Cj/Cq indicate that U-j (r) is repulsive. As can be seen, the achiral 
cis-configuration x = 7r is anti-symmetric with respect to 9 = tt/2 (dash- 
dotted line). This is the case for all ip, and it is a consequence of the additional 
reflection symmetry for a planar configuration. For the other achiral situation, 
X = (the trans-configuration), the potential strength Cj vanishes identically, 
C-j = 0. For all pairs of left- and right-handed configurations, x 7^ 0, 7r, the 
surfaces of equality cannot be transformed into each other by a proper rotation. 
Thus, unlike Uq, the Uj interaction does distinguish between right- and left- 
handed form, the difference depending on both the distance r = £r eq , and the 
orientation f. 
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C 7 (f)/C 6 (f) 




Fig. 4. The relative strength of the discriminative potential C-j (f) compared to the 
dipole-dipole potential Cq (r) in Bohr radii as a function of the polar angle 9 (see 
text). The azimuthal angle is chosen to be ip = ^ and the values of the three different 
twist angles x are indicated in the figure. 



7 Conclusions 



We described a minimalist model for chiral molecules, which, in spite of its 
simplicity, admits a physically plausible and consistent description of their 
chiral properties. As a big advantage of this twisted arc model, the evaluation 
of higher order electric moments can be reduced to a small number of simple 
functions given in closed form. This way the functional dependence of the mul- 
tipole polarizabilities on the frequency and on the model parameters shows 
up transparently. The derived chiral properties, such as the rotatory strength 
and the dispersive interaction potentials, thus display a simple, while physi- 
cally consistent dependence on the model parameters. This was demonstrated, 
specifically for the dihedral angle, by evaluating the chirality-distinguishing 
part of the dispersive interaction with a polarizable atom for molecular pa- 
rameters adopted to the dihydrogen disulfide molecule. The comparison with 
the dominant bulk interaction thus permits to assess to what degree right- 
and left-handed, oriented molecules can be distinguished by the dispersion 
interaction. 

This work was supported by the Emmy Noether program of the DFG. 



20 



References 

[I] I. Tinoco and R. Woody, J. Chem. Phys. 40, 160 (1964). 
[2] H. J. Nolte and V. Buss, Tetrahedron 31, 719 (1975). 

[3] N. L. Balazs, T. R. Brocki, and I. Tobias, Chem. Phys. 13, 141 (1976). 

[4] E. Leuliette-Devin, R. Locqueneux, and J. Tillieu, J. Chem. Phys. 75, 1239 
(1981). 

[5] R. K. Kondru, S. Lim, P. Wipf, and D. N. Beratan, Chirality 9, 469 (1997). 

[6] E. A. Power and T. Thirunamachandran, Chem. Phys. Lett. 3, 361 (1969). 

[7] W. Rijks and P. E. S. Wormer, J. Chem. Phys. 90, 6507 (1989). 

[8] J. A. Cina and R. A. Harris, Science 267, 832 (1995). 

[9] Y. N. Chiu, A. V. Kenney, and S. H. Brown, J. Chem. Phys. 73, 1422 (1980). 

[10] A. J. Stone, The Theory of Intermolecular Forces, Clarendon Press, Oxford, 
1996. 

[II] C. Mavroyannis and M. J. Stephen, Mol. Phys. 5, 629 (1962). 

[12] J. K. Jenkins, A. Salam, and T. Thirunamachandran, Phys. Rev. A 50, 4767 
(1994). 

[13] D. P. Craig and T. Thirunamachandran, Theor. Chem. Acc. 102, 112 (1999). 

[14] E. A. Power and T. Thirunamachandran, Proc. Roy. Soc. Series A 313, 403 
(1969). 

[15] D. N. Sears, C. J. Jameson, and R. A. Harris, J. Chem. Phys. 119, 2694 (2003). 
[16] A. Rauk, J. Am. Chem. Soc. 106, 6517 (1984). 

[17] D. P. Craig and T. Thirunamachandran, Molecular Quantum Electrodynamics, 
Academic Press, London, 1984. 

[18] A. D. Buckingham and M. B. Dunn, J. Chem. Soc. (A) , 1988 (1971). 

[19] M. Pericou-Cayere, M. Rerat, and A. Dargelos, Chem. Phys. 226, 297 (1998). 

[20] A. D. Buckingham, Permanent and Induced Molecular Moments and Long- 
Range Intermolecular Forces, in Intermolecular Forces, edited by J. O. 
Hirschfelder, John Wiley and Sons, New York, 1967. 

[21] L. D. Landau and E. M. Lifschitz, Quantenmechanik, Akademie Verlag, Berlin, 
1988. 

[22] Y. M. Chan and A. Dalgarno, Proc. Phys. Soc. 86, 777 (1965). 

[23] V. P. Osinga, S. J. A. van Gisbergen, J. G. Snijders, and E. J. Baerends, J. 
Chem. Phys. 106, 5091 (1997). 



21 



